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Abstract 

The configuration manifold M of a mechanical system consisting of two uncon¬ 
strained rigid bodies in R", n > 1, is a manifold with boundary (typically with 
singularities.) A complete description of the system requires boundary conditions 
that specify how orbits should be continued after collisions. A boundary condition 
is the assignment of a collision map at each tangent space on the boundary of M 
that gives the post-collision state of the system as a function of the pre-collision 
state. Our main result is a complete description of the space of linear collision 
maps satisfying energy and (linear and angular) momentum conservation, time 
reversibility, and the natural requirement that impulse forces only act at the point 
of contact of the colliding bodies. These assumptions can be stated in geometric 
language by making explicit a family of vector subbundles of the tangent bundle to 
the boundary of M: the diagonal, non-slipping, and impulse subbundles. Collision 
maps at a boundary configuration are shown to be the isometric involutions that 
restrict to the identity on the non-slipping subspace. The space of such maps is 
naturally identified with the union of Grassmannians of fe-dimensional subspaces of 
]^n-i, Q < /j < each subspace specifying the directions of contact roughness. We 

then consider non-standard billiard systems, defined by fixing the position of one 
of the bodies and allowing boundary conditions different from specular reflection. 
We also make a few observations of a dynamical nature for simple examples of 
non-standard billiards and provide a sufficient condition for the billiard map on the 
space of boundary states to preserve the canonical (Liouville) measure on constant 
energy hypersurfaces. 


1 Introduction 

The classical theory of collisions of rigid bodies provides a very natural setting in which 
to explore the geometry and dynamics of mechanical systems on configuration manifolds 
with boundary. From this geometric perspective, the response of the system to collisions 
between its rigid moving parts is specified by assigning appropriate boundary conditions 
that tell how a trajectory should be continued once it reaches the boundary. For example, 
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in the theory of billiard dynamical systems, a topic that may be defined very broadly 
as the study of Hamiltonian (more typically, geodesic flow) systems on Riemannian 
manifolds with boundary, one typically assumes that trajectories reflect off the boundary 
specularly—the simplest form of impact response compatible with the basic laws of 
mechanics such as energy conservation and time reversibility. Billiard systems with more 
general boundary conditions have to our knowledge been investigated very rarely. One 
pertinent example is |3], which is restricted to 2-dimensional billiards. (There is, of course, 
an extensive literature in engineering and applied physics about less idealized systems 
governed by impact interactions, but this literature is not concerned with the differential 
geometric issues that are the main focus here.) 

Our first goal is to classify boundary conditions for systems defined by two unconstrained 
rigid bodies in K", n > 1, under standard physical assumptions of energy conservation, 
linear and angular momentum conservation, time reversibility, linearity of response, and 
another condition to be defined shortly that extends momentum conservation and is 
typically made implicitly in textbooks. Collisions satisfying all of these properties will be 
called strict. They are formally represented by linear maps Cq : TqM TqM, where M 
is the configuration manifold equipped with the kinetic energy Riemannian metric, q is 
a boundary configuration, and the tangent space TqM is the space of (pre- and post-) 
collision states. A boundary condition for the system then consists of the (differentiable, 
measurable, random, etc.) assignment of a strict collision map Cq to each boundary 
configuration q e dM. 

In dimensions greater than 1, the collision map is not uniquely determined by the 
conditions of strict collision. It is well-known that the nature of the contact between the 
colliding rigid bodies also needs to be specified. The standard case in which Cq is specular 
reflection corresponds to bodies having physically smooth surfaces. 

Towards this classification we identify a family of subbundles of T{dM) arising naturally 
under the assumed physical laws and discuss some relationships among them. We then 
show examples of trajectories of non-standard (i.e., non-specular) boundary conditions and 
make a few observations about their dynamics based on numerical simulation. Although 
we leave for future work a more systematic analysis of rough billiard dynamics, we give 
here sufficient conditions for the non-standard billiard system to leave invariant the natural 
volume measure on a constant energy manifold (derived from the canonical symplectic 
form). The invariance of this billiard measure makes it possible to bring the tools of 
ergodic theory (see [11|8]) to the study of non-standard billiard systems, although we do 
not pursue this direction here. 


2 Statements of the main results and examples 

2.1 Notation, terminology, and standing assumptions 

For the most part we consider the unconstrained motion of two rigid bodies, represented 
by the sets Bi,B 2 c K". We call these sets the bodies in reference configuration. Let 
G = SE{n) denote the Euclidean group of orientation preserving isometries of R" equipped 
with the standard inner product. The bodies are assumed to be connected n-dimensional 
submanifolds of R" with smooth boundaries. An interior configuration of the system 
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consists of a pair (51,(72) ^ G x G such that 5i(i?i) and 52(^2) are disjoint sets. The 
closure of the set of interior configurations, denoted M, will be called the configuration 
space of the system, and its boundary is the set of contact (or collision) configurations. 
M has dimension 2 dim(G) and the nature of the boundary dM will depend on geometric 
assumptions about the Bj. We will soon state a sufficient, and fairly general for our 
needs, condition on the dBj for M to be a submanifold of G x G with smooth boundary. 

The following notations are used fairly consistently throughout the paper. Points in 
Bj are denoted b,bj. Elements of the Euclidean group, which is the semidirect product 
G = SO{n) X M” of the groups of rotations and translations, are written as pairs {A, a), 
where A 6 SO{n) and a 6 M". Elements of the Lie algebra g = se(n) are written {Z,z), 
possibly with subscripts or superscripts, where Z 6 so{n) and z 6 M". The outward¬ 
pointing unit normal vector to the boundary of Bj at b € dBj is denoted Vj{b). It is 
convenient to consider orthonormal frames a at b e dBj adapted to the bodies, in the 
following sense: cr: M" ^ TbM" = M" is an element of SO{n) such that CTe„ = -(- 1 )^i'j{b), 
where e„ = (0 ,..., 0,1)^ is the last element of the standard basis {ei,..., e„} of M" and 
‘f ’ indicates matrix transpose. 



Figure 1 : On the left and right are the bodies Bi and B2 in their reference configuration in R". 

A configuration of the system of rigid bodies is given by a pair (51,(72) of elements of 
the Euclidean group G. A boundary configuration can be parametrized by the tuple 
(61,(71,7)2,(72) where bj e dBj such that 51(61) = 52(62) and aj is an orthonormal 
frame at bj as will be explained in the text. 


Because M is a submanifold of G x G, each tangent space can be canonically identified 
with T^e,e) {G X G) = g ® 0 by left-translation. Thus states of the system, defined as 
elements of TM, may be canonically identified with tuples (Ai, ai, A2,02, .^1, ^i, .^2, ^2)- 
We sometimes indicate the state by ((7,^), where q = (51,52) and ^ = {Zi,Zi,Z2, Z2) 6 0X0- 
The position of material point b e Bj in the given state is then Pjfib) = Ajb + aj and 
its velocity is V{b) = Aj{Zjb + Zj). The boundary configuration (51,52) can also be 
parametrized, up to an overall rigid motion of the two bodies keeping their positions 
relative to each other unchanged, by (&i,(ii,62,(72), where bj is in the boundary of Bj 
and aj is an adapted frame such that 51(61) = 52(62) and Aiai = A2CT2. The tuple 
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(6i,(Ti/i,62,(72^) corresponds to the same contact configuration, for all he H, where 


H := SO{n - 1 ) = {A 6 SO{n) : = e„}. 

These notions are illustrated in Figure 

Let Ilh denote the orthogonal projection from M” to the tangent space to the boundary 
of Bj at a boundary point b. The orthogonal projection to will be denoted 11 . 

The shape operator of the boundary of Bj at the point b is the linear map defined by 

* F? ~ ~ByVj 

where Dy is directional derivative in K". We say that S ■ ^ is the shape 

operator Sb in the adapted frame cr at 6 if aS = Sbcr- The notation Ado-iS" = aSa^^ will be 
used often to indicate conjugation. 

So as not to get distracted by regularity issues, we assume that the configuration 
manifold M has smooth boundary and that each boundary configuration corresponds to 
the bodies being in contact at a single common point. Proposition | 2 . 1 [ which is a special 
case of Proposition | 5 .H gives a sufficient condition for M to be nice in this respect. The 
hypotheses of the proposition will be assumed to hold throughout the paper. 

Proposition 2.1. Suppose that Si + S2 is nonsingular for every relative configuration 
( 5 i,(Ti,62,(72) of the rigid bodies, where Sj : -»• is the shape operator of 

the boundary of Bj at bj in the adapted frame Uj. Then M is a smooth manifold of 
dimension 2 dimG' with smooth boundary, and each boundary point q = (51,52) represents 
a configuration with a unique point of contact. Moreover the map that associates to 
q 6 dM the uniquely determined pair (61,62) € dBi x dB2 such that 51(61) = 52(62) is 
smooth. 


We call the 61,62 associated to 5 € dM under the condition of Proposition 2.1 the 
contact points (in the reference configuration) associated to boundary point q. 


2.2 The kinematic bundles 

If a, 6 6 K", let a a 6 € so{n) be the n-by-n matrix such that (a a b)ij = Ujbi - aibj. If a, 6 
are orthogonal unit vectors, a A 6 is the infinitesimal generator of the one-parameter group 
in SO{n) that rotates the plane spanned by a and 6 and fixes pointwise the orthogonal 
complement of that plane. The boundary state of the two-body system consists of the 
boundary configuration q = (51,52) e G x G and velocities ^ = {Zi, zi, Z2, zi) e g x g. 

We now define the kinematic bundles. Given q = (51,52) e dM, with gj = {Aj,aj) and 
associated contact points 61,62, consider the following linear relations on the ^ € Tq{dM), 
where Nj = dBj and Vj = i'j{bj')\ 

Ri : vi ■ {Zibi + zi) = V 2 • (.^262 + Z2) 

R 2 ■ Ai{Zibi + Zi) = A2{Z2b2 + Z 2 ) 

i ?3 : Ad^^ Zj = W + Vj A Wj for W 6 so(n) and Wj e Tb^ A^, j = 1,2 
R 4 '. AdiA^Zi = Ad^2-^2- 
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As already noted, Aj{Zjbj + Zj) is the velocity of the contact point bj in the given state. 
Observe that relation R2 implies Ri. It will be shown later that 

Tg{dM) = {CeQXQ:R^}. 

The physical interpretation of the kinematic bundles is as follows. A state satisfying 
i?i has the property that the contact points have zero relative velocity in the normal 
direction to the plane of contact. Relation R2 is satisfied exactly when the contact points 
are not moving at all relative to each other at the moment of contact. This means that 
the contact points do not slip past each other. Relation R3 describes a state in which the 
tangent spaces to the bodies at the point of contact do not experience a relative rotation 
(on that tangent space). Thus it is a condition of non-twisting. Together R2 and R^ 
describe a state in which the bodies are rolling on each other. 

Definition 2.1 (Kinematic bundles). Let &, 91 , and © be the vector subbundles of 
T{dM) defined by 


&q = {f€QXQ-.R2} 

91 , = {Ce 0 xg:R 2 ,R 3 } 

Dg = {f€QXQ:R2,R3,R4}. 

Note that 29 c 91 c ©. We refer to © as the non-slipping subbundle, 91 the rolling 
subbundle, and 29 the diagonal subbundle. 

It will be shown that the diagonal subbundle 29 is the tangent bundle to the orbits of the 
action of G on M by left translations: 5(51,52) •= (5511552)- Later in we give a different 
definition of these subbundles. Definition | 5 . 2 [ that makes their physical interpretation 
more clear. Then what is stated above as a definition is derived in Section [531 


2.3 The kinetic energy metric and the impulse subbundle 


Suppose now that the bodies Bj are assigned mass distributions represented by finite 
positive measures g.j supported on Bj. Let ruj := fj,j{Bj) be the mass of Bj. We may 
assume without loss of generality that g,j has zero first moment: fg bdfj.j{b) = 0 . This is 
to say that Bj has center of mass at the origin of K". The matrix of second moments of 
fj,j is Lj = (Irs), with entries 


I 


rs 



We call Lj the inertia matrix of body Bj. This matrix induces a map Lj on so(n) that 
associates to Z e so(n) the matrix Lj{Z) = LjZ + ZLj € so(n). 

Definition 2.2 (Kinetic energy Riemannian metric). Given q e M and u,v e TqM, define 
the symmetric non-negative form on TqM by 


{u,v)q = Y,mj 
3 


TT(L,{Z^)Z]f) 


+ z. 
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where 2;“, ^2 , Z2 ) (^1 j,-^2 j- 22) are the translates io g x g ofu,v. When the 

above bilinear form is positive definite we call it the kinetic energy Riemannian metric 
on M. Denoting by || • ||q the corresponding norm at q, we call |||'y||q the kinetic energy 
associate to state {q,v). 


The kinetic energy function given in Definition 2.2 is easily shown (as indicated later) 
to come from integration with respect to the mass distribution measures of (one-half of) 
the Euclidean square norm of the velocity of material point b over the disjoint union of 
Bi and B2■ Thus Definition | 2 . 2 | agrees with the standard textbook definition of kinetic 
energy. It is also clear that the metric is invariant under the left-action of G on M. Note 
that the boundary of M is a C?-invariant set. 

For each it 6 g we define vector field q*^ UqC TqM, q = (51,52), by 


We call u >-->■ u the infinitesimal action derived from the left G-action on M and u the 
vector field associated to it 6 g. 

Definition 2.3 (Momentum map). The map IP® : TM g* defined by 


T®(5,5)(m) = {q,Uq)q 


is called the momentum map associated to the G-action on M. 

The most straightforward way of introducing dynamics into the system is through 
Newton’s second law. There are several equivalent forms of it as we note later. The 
following is particularly convenient for our needs. We first define a force field (possibly 
time dependent) as a bundle map F : TM T*M. Given a state (5,5), q = (51,52), each 
component Fj of F can be pulled-back to g* using right-translation Rg^, so it makes sense 
to write 

( 2 . 1 ) 

This is Newton’s second law written as a differential equation on the co-Lie algebra of G. 
Other useful forms are mentioned later. One of them is indicated in the next proposition, 
in which we use the notation F"^ for the dual of F with respect to the left-invariant 
Riemannian metric and write 


{Yj(t,q,q),yg{fiq,q)) := {dLg^)^^F*{t,q,q) e q. 

Here we are using the differential of the left-translation map Lg.. 

Proposition 2.2. The equation ^J’j( 9 , 5 ) ~ YGg-Fj is equivalent to 

mj {dLjZj - [LjZj, Zj^'j = Lj{Yj) 
mjVc = AjPj 

where gj = and Vc = AjZj is the velocity of the center of mass of body Bj. 
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We assume that F results from the integrated effect of forces acting on the individual 
material points. That is, we assume that there exists a IR"-valued measure ipj on Bj 
parametrized by TM from which F is obtained by integration: 

F{q,v){u)= f Vu{b) ■ d(pj^q^y{b) 

J Bj 

for all u € TqM, where Vu{b) is the velocity of the material point b in the state {q,u). Of 
special interest for us are the forces involved in the collision process. These impulsive 
forces are characterized by being very intense and of very short duration, applied on a 
single point—the point of contact in each body. 

That the forces act on each body only at the point of contact greatly restricts the 
right-hand side of the equation of motion in Proposition | 2 . 2 [ This is indicated in the next 
proposition. 

Proposition 2.3. We suppose that the force field Fj acting on body Bj is such that the 
force distribution measure ip is singular, concentrated at the point bj. Then the equations 
of motion of Proposition \ 2 ^ reduce to 

mj {^LjZj - ^LjZj, Zj\^ = bj a pj 

mjVc = AjPj 

For ideal impulsive forces (of infinite intensity and infinitesimal duration), momentum 
should change discontinuously. Integrating Equation | 2 . 1 | over a very short time interval 
[F,F] around t produces a nearly discontinuous change in momentum while keeping the 
configuration essentially unchanged. We have informally 

3 ^® (9, g+) - J’® (9,) = ^ Fj ds = Impulse at t. 

It is not necessary for our needs to make more precise the limit process suggested by this 
expression. From it we obtain the form of the change in momentum after impact, which 
is given in the next proposition. Let q = (31,52) £ dM be a collision configuration and 
denote by 

{Zf,zt,Z^,z^)eTqM 

the post- (+) and pre- (-) collision velocities of the two rigid bodies. 

Proposition 2.4 (Velocity change due to impulse at contact point). Given pre-collision 
velocity {Zf,zf,Z2,zf^) there exist ui,U2 e K" such that 

^j = ^j + 

Zj = Zj + Tj^{bj A Uj). 

Under conservation of linear momentum miAiUi + m2A2U2 = 0 holds. 

The proof of the above proposition is given in Section | 4 . 4 | The assumption that 
impulsive forces of one body on the other at the moment of impact are applied at the 
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point of contact is a strong constraint. One can in principle conceive of force fields of 
relatively long range, acting throughout the bodies, that are briefly switched on at the 
moment of impact, then switched off as soon as the bodies lose contact. More realistically, 
the bodies could suffer a deformation around the region of impact, creating a small 
neighborhood of contact. Of course this goes beyond the rigid body model. Here it is 
assumed that these possibilities do not happen, and that any effect of one body on the 
other can only be transmitted through the single point of contact between them. 

If Lj is non-negative definite of rank at least n- 1 , Lj is invertible. With this in mind, 
Proposition | 2 . 4 | suggests the following definition. 

Definition 2.4 (Impulse subbundle). The impulse subbundle ofTM (over the base 
manifold dM) is defined so that its fiber at q € dM is the subspace 

^ ui), ui), (£ 2 ^{b2 A M 2 ), M 2 )) : Uj 6 ,miAiUi + 1712A 2 U 2 = 0 } . 

We have now the following vector subbundles of i*{TM), where i : dM M is the 
inclusion map: S) c 91 c 6 c T{dM) and £. The latter subbundle is the only one that 
depends on the mass distributions. 

Theorem 2 . 1 . The impulse subspace €q is the orthogonal complement of the non-slipping 
subspace <Sq and contains the unit normal vector n^. Therefore, 

TqM = &q® (€q Q Kllg) ffi 


is an orthogonal direct sum. 

A physical interpretation of this orthogonal decomposition will be given shortly. 

2.4 Collision maps 

Let Hq be the unit normal vector to dM pointing into M at a boundary configuration q. 
Define the half-spaces 


T;M := {v 6 TqM : (M,n,) > 0 } = -T'M. 

We call any Cg : T“ M ^ TqM a collision map at q. By a boundary condition we mean 
the assignment of such a map C, to each q € dM. We only consider here linear collision 
maps; that is, Qq extends to a linear map on TqM. 

Definition 2.5 (Strict collision maps). A collision map Cq at q e dM is strict if the 
following hold for all u,v € TqM: 

1 . Conservation of energy.- {CqV,Gqu)q = {v,u)q. That is, Gq is a linear isometry. 

2 . Conservation of momentum; ‘J’^{q,Gqv) = y^{q,v). 

3 . Time reversibility; = Id (a linear involution). 

4 -. Impulse at the point of contact; GqV - v e Gq. 


Proposition 2.5. Condition 2 of Definition \2. ^ is equivalent to assuming that Cg restricts 
to the identity map on iDq. Condition f is equivalent to CgCg = Cg and (Cg - Id)£q = 0. 
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Thus energy conservation and impulse at a single contact point are together equivalent 

g. In this sense, condition (4) of Definition 


q being the identity on &q 


2.5 


can be 


regarded as generalizing momentum conservation as we irote in Proposition |2.5| In fact, 
coirservatioir of momeirtum amounts to Gg beiirg the ideirtity oir T>g, whereas 4 aird 
Theorem 2.1 imply that Gg is the identity on the bigger subspace &g. An intermediate 
condition is that Gg restricts to the identity on the rolling subspace 9Ig. 


Corollary 2.1. Strict collision maps are the linear isometric involutions ofTgM, q € dM, 
that restrict to the identity map on the non-slipping subspace &g. 

Collision maps have eigenvalues ±1. The map := is the orthogonal projectioir 
to the eigenspace associated to eigenvalue ±1. 


Definition 2.6. The dimension of the eigenspace of Gg associated to eigenvalue -1 will 
he called the roughness rank of Gg. The image of the orthogonal projection P- (a subspace 
of Gg) will be called the roughness subspace at q. 

The unit normal vector is always contained in the impulse subspace Gg and it must 
irecessarily be in the - 1-eigenspace of Gg. 

Corollary 2.2. Let n he the dimension of the ambient Euclidean space. Identifying 

Gg © MlUq = ^, 

the set of strict collision maps is the set of C e 0{n- 1) such that = I. Writing 


3k ■■= 0{n-l)/{0{n-k-l)xO{k)), 


then the set of strict collision maps at any given boundary point is 3o^' "^3n-i- Moreover, 
dimSk = k{n- k-1) and k is the roughness rank at q. We call 3k the Grassmannian of 
rough subspaces having roughness rank k. 


It is easy to compute the dimensions of the Grassmannians 3k for strict collision maps. 
They are giveir, up to dimension 5, by the followiirg table: 


k 

n 

0 12 3 4 

1 

0 

2 

0 0 

3 

0 1 0 

4 

0 2 2 0 

5 

0 3 4 3 0 


The table shows that in dimension 1 there is a unique strict collision map; in dimension 2 
there are exactly 2 possibilities; and in dimension 3 there is one possibility of roughness 
rank 0 given by the standard reflection map, one possibility for maximal roughness rank 2, 
and a one-dimensional set of possibilities for roughness rank 1 parametrized by the lines 
through the origin in For geireral n, the unique collision map of maximal roughness 
rank will be referred to as the completely rough reflection map. 
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2.5 Non-standard billiard systems 


We have so far considered systems consisting of two unconstrained rigid bodies. The 
results of this paper can be extended to situations in which one body or both are 
subject to holonomic and non-holonomic constraints. We will explore this extension more 
systematically elsewhere. Here we consider only the case in which body Bi remains fixed 
in place whereas B 2 is unconstrained except for the condition that it cannot overlap with 
Bi. The term billiard system will refer to a system of this kind where B 2 is a ball with 
rotationally symmetric mass distribution. The system will be called non-standard if the 
(strict) collision maps are not all specular reflection. 


reference configuration 



Figure 2 : A billiard system. Body B\ is kept fixed (the billiard table) and B := B2 is a ball of 
rotationally symmetric mass distribution that can move freely in the complement of 
Bi. Given the contact point b s N and the post-collision velocity v of the center of 
mass of B, the point of contact of the next collision will be written b' = T{b,v). 


Let R denote the radius of B := B 2 and m its mass. Due to rotational symmetry of the 
mass distribution of B, the matrix of inertia L is scalar, that is, L = A/. For example, a 
simple integral calculation shows that if B has uniform mass distribution, then A = . 

(Recall that we have defined L as the matrix of second moments of the mass distribution 
measure divided by the total mass, so m does not appear in A.) The (smooth) boundary 
of Bi will be denoted N and the unit normal vector field on N pointing into the region 
of free motion of B will be denoted v. Trajectories of the billiard system are sequences of 
states: ( 50 , ?o), ( 51 , Ci), • • ■, where 

( 91 ,^ 2 ) ^SE{n)xst{n) = TSE{n), gi = {A„ai), = 

Here gi is the contact configuration and the post-collision velocities in the body frame 
(reference configuration) at the ith collision. 
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To each contact state {g,0 = (^jO, Z, 2 ) is associated a unique contact point h e N 
and the post-collision velocity v = Az of the center of mass. The center of mass of B in 
configuration g is a and the velocity of any given material point b € B is V (5) := A{Zb + z). 
When it is necessary to distinguish points in N and in B we write b e N and ba € B. The 
unit normal vector to B at bo will be written lyo(bo) = bojR. The point of contact at the 
next collision, which only depends on b and v, will be denoted b' = T{b,v). See Figure]^ 

One step of the billiard motion, (g,^) 1 -^ {g',^'), amounts to the following operations. 

1. From the current collision state (g,^) at time t one obtains the contact point & € TV 
and velocity v of the center of mass a oi B where g = {A, a) = {A{t),a{t)). It should 
be kept in mind that ^ = {Z,z) describe post-collision velocities so v = Az points 
into the region of free motion of the ball. 

2. Obtain the contact point b' = T{b,n) & N and the time t' = t + r ot the next collision. 

3. Obtain the next pre-collision state: {g',^~) where g' = {A',a') = (A(t + r), a(t + r)), 

= (Z~,z~), and 

A' = Ae^^, a = a + tAz, Z~ = Z, z~ = e~'^^z. 

This is the free (geodesic) motion between collisions. Observe that a' = b' + Rv{b'). 

4. Let bo = {g')~^b' e dB be the contact point on the ball in the reference configuration 
at the next collision and denote by Ho,!!;) the orthogonal projections to the tangent 
space to dB at bo and to 'Ri'o{bo), respectively. Note that A'vo{bo) = -v{b'). 

5. Finally, compute = {Z'^z') from {Z~^z~) using the choice of collision map. It 
will be shown that 

(2.2) (Z', z') = [z--^boA{I- 7)V-,z- - ail - 7)^“ - 2niz-) , 

where a := 1/(1 + R^I2X), V~ = Ilo(Z~bo + z~), and T is a linear involution on 
(dB) corresponding to a choice of collision map. For specular reflection 7=1 
and for completely rough collisions T = -I. 

2.6 Examples of non-standard billiards 

We assume in all examples the uniform mass distribution on the ball B so X = R?Hn + 2), 
where R is the radius of B. Let first n= 2. In this case the only non-standard collision 
map corresponds to T = -I. Elements of the rotation group are parametrized by the angle 
of rotation 0 and elements of the Lie algebra of iS'0(2) are written as 0J, where J is the 
rotation matrix by 7 r/2 in the counterclockwise direction. Together with the standard 
coordinates ix,y) we obtain coordinates {9,x,y) on SE{2). It will be convenient to 
make the coordinate change: xq = ROj\/2,xi = x,X 2 = y- This yields coordinates 
(xq, xi, X 2 , xq, xi,X 2 ) on the billiard state space. We also write vq = xq and v = (ii,i 2 )^ 
for the velocity of the center of mass of the disc. 

The choice of coordinates is made so that the kinetic energy Riemannian metric becomes, 
up to multiplicative constant, the standard Euclidean metric. Then it can be derived 
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from Equation 2.2 that the post-collision velocities (vg, 
contact be N is the function of the pre-collision velocities 


v^) after collision at point of 
(^' 0 ;^^”) given by 


(2.3) 


^ 1 _ 2^2 , ,,,, 
^0 =“3^0 + 

V = -^1^0 + 


Jv{b)-v ■v{b)v{b). 


Thus the state updating equations for a 2-dimensional non-standard billiard system is as 
follows. If T is the time of free flight between the two consecutive collisions and setting 
X = {xg,Xi,X 2 ), V = (uo,Ui,W 2 )) then the billiard map giving the next state (x', v') as a 
function of the present state (x,v“) is (x',v') = (x+ rv,v’*') where is related to v“ 
according to Equations |2.3| The geometric interpretation of those equations is explained 
in Eigure|^ Note the role played by the angle /3 defined by cos/3 = 1/3, sin/3 = 2\/2/3. 
In |3] it is observed that /3 is the dihedral angle of a regular tetrahedron. In the figures 
to follow we only indicate the position of the center of the disc; we draw a smaller table 
whose boundary is at a distance R from the boundary of the original table and we imagine 
the center of the ball as a point mass bouncing off the boundary of this smaller region. 





-Ju{b) 


Figure 3: Geometric interpretation of the rough reflection in angle-position space for n = 2. The 
eo component of the incoming velocity v is the scaled angular velocity vo = i?@/\/2. 
The outgoing velocity is obtained by first reflecting v specularly on the plane spanned 
by eo and Jv(v) to find v, then reflecting the latter specularly on the plane spanned 
by v{b) and Ju(b), and finally rotating the resulting vector by /3 as indicated. As 
noted in [^, d is the dihedral angle of a regular tetrahedron. 


Next we show examples of trajectories of systems with rough collisions. The examples 
are given here without much analysis. We leave the more systematic study of the dynamics 
of such systems for another article. As a first illustration, consider the case of a circular 
billiard table. The typical trajectory is shown in Figure]^ and a few more examples are 
shown in Figure 

The following proposition captures the main properties of trajectories of circular rough 
billiards readily observed in the Figure]^ 
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Figure 4: Caustics of a circle billiard with rough collisions consist of pairs of concentric circles. 

The angle a at each vertex of the projection of a trajectory on the xj/-plane is constant 
along the trajectory. 


Proposition 2.6 (Circular billiard with rough collisions). For a billiard system with 
circular table of radius r and rough collisions, the projections of trajectories from the 
3-dimensional angle-position space to the disc in position plane have the property that 
the vertex angle at each collision is a constant of motion. Moreover, for each projected 
trajectory 7 , there exists a pair of concentric circles of radius less than r that are touched 
tangentially and alternately by the sequence of line segments of j at the middle point of 
these segments. 



Figure 5: Three orbit segments with different initial conditions for the motion of the center of 
mass of a disc in a circular billiard table with rough contact. 


The next example consists of a moving disc in a wedge-shaped table with rough collisions. 
A few examples of trajectories for different values of the vertex angle of the billiard table 
are shown in Figure]^ What is most notable in this case is the existence of bounded 
orbits. Other properties such as periodic orbits for certain angles of the wedge table and 
caustics are clearly suggested by the figures. 
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Figure 6: Motion of a disc in wedge shaped table. The typical segment of trajectory (more 
precisely, the motion of the center of mass only) is shown on the right. The two 
trajectories on the left, which consist of 1000 free flight segments each, are likely 
periodic. 


In the previous examples the boundary condition on M amounted to a constant (more 
precise, parallel) choice of Qq. We wish to illustrate now boundary conditions for which 
the map q >->■ Gq varies in a nontrivial way or is chosen randomly. Let the billiard system 
consist of a disc moving in an infinite strip bounded by two parallel lines. We suppose 
that one hemisphere of the boundary of the disc is rough and the other is smooth. In 
Figure]^ we show graphs of the position of the (center of) the disc along the longitudinal 
axis of the table as a function of the collision step. The time between two consecutive 
collisions is easily shown to be constant, so the step number is proportional to time. The 
three graphs describe the same trajectory at different time scales, as indicated in the 
legend of the figure. 



Figure 7: A single orbit of the motion of a disc between parallel plates in dimension 2. Half 
of the boundary circle is rough and the other half is smooth. The horizontal axis 
indicates the step number, taken as a proxy for time. The vertical axis gives the 
distance of the center of mass along the length of the 2-dimensional channel. 


It is interesting to observe the apparent long range quasi-periodic behavior of trajectories. 
It is also interesting to note the differences between this example and the next shown in 
Figure The setting is essentially the same, except that a point on the boundary of the 
disc is chosen to be rough or smooth randomly with equal probability. This is thus an 
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example of a random boundary condition. The longitudinal motion now corresponds to a 
random walk, for which it is possible to prove a diffusion (Brownian motion) limit. 



Figure 8: Disc between parallel plates in dimension 2. This is similar to Figure [^except that 
the roughness rank is now random, either 0 or 1 with equal probabilities. 


We consider now a few examples in dimension 3. In all cases, a ball of uniform mass 
distribution moves between two parallel infinite plates. In dimension three, the roughness 
rank can be 0, I, or 2. Standard specular reflection has roughness rank 0; we explore 
examples of roughness rank 1 and 2. 



Figure 9: Typical segment of trajectory in position space and its horizontal projection for the 
motion of the center of mass of a ball bouncing between two parallel plates in dimension 
3 with roughness rank equal to 2. Note that orbits are bounded. 


Figure illustrates the case of roughness rank 2 collisions. The figure on the right 
shows the projection of the trajectory to the coordinate plane parallel to the plates. One 
notable property of the system is that trajectories are bounded. In dimension 2, a similar 
property was noted in |3]. 

When the roughness rank is one, the set of collision maps comprise a one-dimensional 
family. Figure^] shows some combinations of boundary conditions. One clearly notices 
that whenever the roughness rank in at least one plate is not maximal, trajectories are 
no longer bounded. The boundary conditions for the systems of Figure [T0| are of the 
following types: one plate has roughness rank 2 and the other has roughness rank 1 with 
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constant rough direction (that is, constant map T in Equation 2.21; and both plates have 
roughness rank 1, with random rough direction for the bottom plate and either constant 
or independent random rough direction for the top plate. Specifically, we choose the 
directions given by angles 0,7r/3,27r/3 with equal probabilities. The legend of the figure 
shows which trajectory corresponds to which condition. 




Figure 10: Horizontal projection of motion of the center of mass for a ball bouncing between two 
parallel plates in dimension 3. Top left: roughness rank 1 for both top and bottom 
plates and random (and independent) roughness directions; bottom left: top and 
bottom roughness rank 1, but now roughness direction is constant for the top and 
random for the bottom plate; right: roughness rank is 2 for top plate and 1 for bottom 
plate, with constant roughness direction. Making the the rough direction random for 
the bottom plate gives a trajectory that does not look significantly different than 
the one on the top left. 


2.7 INVARIANT MEASURE 

A fundamental property of the dynamics of standard billiard systems is the existence 
of a canonical invariant measure on constant energy surfaces, sometimes referred to as 
the Liouville measure. We give here a sufficient condition for the same measure to be 
invariant under non-standard collisions. 

Let S denote the boundary of the configuration manifold of the two-bodies system. As 
before, we assume that S is smooth. We fix a value £ of the kinetic energy and denote 

N^ = i^{q,v)€TM:qeS,^\\vf = £.y 

Define the contact form 9 on TM to be the 1-form such that 9y{^) = {v,dTTy^)q, where tt 
is the base point projection from TM to M (and we indicate the element of TM by u 
rather than {q,v) in subscripts). It is well-known that dO defines a symplectic form on 
TM. It can also be shown that the restriction of d9 to defines a symplectic form on 
\ TS. (See, for example, |S].) 
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The billiard map T on associates the post-collision state of the system at the 
time of a collision to the post-collision state at the next collision. There are well-known 
issues about this map, even for standard billiards in dimension 2, that make the precise 
specification of its domain difficult to describe. See, for example, [1|. Here we assume 
that the domain of T consists of a “large” open set of full Lebesgue in and omit any 
further reference to it since this issue of domains is not specific to our rough billiards. 
The next result is shown by a local argument and considerations of domain do not play a 
role. 



Figure 11: Angle-position parallelepiped for rectangular billiard table and disc. The canonical 
invariant billiard measure on a constant energy hypersurface is, up to multiplicative 
constant, the product of the Euclidean area measure on the boundary of 31 and the 
measure on the hemisphere of velocity directions given by cos (j) dA where dA is the 
Euclidean area measure on the hemisphere. If billiard trajectories are initiated on 
the side X 2 = 0 with random initial condition given by the jus t de scribed measure, 
return trajectories will have the same distribution. See Figure 
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Theorem 2.2. Suppose that the field of collision maps q € S Gq is piecewise smooth 
and parallel (where it is smooth) with respect to the Levi-Civita connection associated to 
the kinetic energy Riemannian metric. Let LI = d9 a ■■■ a dO he the form (of degree 2n - 2, 
where n is the dimension of the ambient Euclidean space) derived from the canonical 
symplectic form dO on \TS. Then LI is, up to sign, invariant under the billiard map. 

Corollary 2.3. Rough billiards in dimension 2 preserve the canonical billiard measure. 

The theorem will be proved in Sectionj^ Corollary |2.3| is due to the following observation. 
The boundary of M is a flat surface with the Euclidean metric and the vectors cq, Jv{b) 
shown on the right-hand side of Figure [^constitute a parallel frame. The orthogonal line 
distributions Gq and &q are also parallel as the angle between each of them and eo is 
constant. But these are the eigenspaces of Gq for the eigenvalues -1 and 1, respectively. 
It follows that the field of rough collision maps is parallel. 

We will leave for a future paper a more detailed investigation of invariant measures 
of non-standard billiards. Here we simply illustrate Theorem |2.2| with a numerical 
observation concerning the motion of a disc in a rectangular table with rough collisions. 
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The geometric set-up is shown in Figure 11 The configuration manifold in this case is a 


parallelepiped IR in dimension 3 and the canonical billiard measure on the manifold N of 
unit length vectors with base points on the boundary of IR has density proportional to 
p{v) = V ■TUq = cos (j), 0 < (/) < 7r/2, with respect to the Riemannian volume measure on N, 
where (/) is the angle the vector v makes with the normal vector to the boundary. 

It can be shown that Theorem 2.2 applies to this case. As an experiment to illustrate 
invariance of the billiard measure for rough collisions we sample initial conditions on the 
face X 2 = 0 with the uniform distribution for the (ccojCCi) positions and initial velocity v 
having probability density proportional to cos 0 relative to the uniform probability on 
the unit hemisphere. If the return states to the face a ;2 = 0 are distributed according to 
the same measure, then the angle 0 for the return velocity must be distributed relative to 
Lebesgue measure on [0,7r/2] with density sm(2(f>), which is the marginal density function 
for the angle distribution with respect to the Lebesgue measure d(j) on [0,7r/2], under the 
assumption that the billiard measure is invariant. 



0.0 0.5 1.0 1.5 


Figure 12: Experiment to illustrate invariance of the billiard measure for a rectangular table. 


This is indeed the case as shown in Figure 12 A large number (10°) of initial conditions 


starting from one side of the rectangle are sampled from the normalized billiard measure 
restricted to that side. For each trajectory, the return state to that side is computed 
and the angle relative to the normal (to the side of the angle-position parallelepiped 
corresponding to that side of the rectangle) is recorded. The distribution of values is 
shown in the above histogram. The superimposed line is the graph of sin(2(/)). 


3 The Euclidean group and its Lie algebra 

The proofs of the main statements made above will be given after we establish some 
basic material. Despite the classical nature of the subject, standard textbook treatments 
of collisions of rigid bodies are not adequate for our needs while the more differential 
geometric texts in mechanics mostly do not treat this topic. Thus we find it necessary 
to develop the subject more or less from scratch. In this section we review general facts 
about the Lie theory and Riemannian geometry of the Euclidean group with left-invariant 
metrics that will be needed thoughout the paper. 
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3.1 Generalities 


The isometry group of (M”,-), where indicates the standard inner product, is the 
Lie group of all the affine maps of the form x i->- Ax + a for A e 0{n) and a e K" 
under composition of maps. The closed subgroup of orientation preserving isometries, in 
which A € SO{n), is the Euclidean group in dimension n, denoted SE{n). The latter is 
isomorphic to the semidirect product SO{n) x M” with multiplication operation 


(A2, a 2 )(^i, fli) - {A2A1, A 2 ai + 02) 


and inverse 

{A,a)-^ = {A-\-A-^a). 

It is also isomorphic to a subgroup of the general linear group GL{n+ 1,IR) under the 
correspondence 

€GL(n+l,M). 

The Lie algebras of SO{n) and SE{n) will be denoted so{n) and se(n). The former 
consists of all the skew-symmetric matrices in the linear space M(n,M) of n x n real 
matrices and se(n), when SE{n) is viewed as a subgroup of GL{n+ consists of the 
matrices 

where X 6so(n) and x is any vector in M". Indicating the matrix by the pair (X,x), the 
Lie bracket is written 


[{X,x),{Y,y)] = {XY -YX,Xy -Yx). 
One-parameter subgroups of SE{n) have the form 


ojj'io 1 J' 


It is useful to introduce the wedge product, the bilinear operation that associates to a 
pair of vectors a, b in K" the skew-symmetric matrix a a b e so(n) whose (z, j)-entry is 

(a A b)ij = Ojbi - Gibj. 

The following elementary properties of the wedge product will be used. The transpose of 
a matrix will be indicated by . 

Proposition 3.1. Let a,b,u be (column) vectors in ISA , AeSO{n) andZeso{n). Then 

1. (a Ab)u = (a ■ u)b - (b ■ u)a 

2. (a Ab)f = b A a 

3. A{a A b)A~^ = (Aa) a (Ab) 
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4- Tr ((a A = 2(Za) • 6 

5. Tr ((a A &)(c A d) Q = (a • c)(&• fi) 

6 . Let V be the span of orthogonal unit vectors a,b e M". Then (a a b)^ = -I and 

R{0) ■= exp(da a &) = (cos0)I + (sind)a a 6 6 SO{n). 

Thus R{0) is the identity on , and a rotation on V. 

7. Let e„ = (0,..., 0,1) ^ € M" and 11: M” -»• = ef the orthogonal projection. Then 

Z = UZn + e„ A (Zcn) 

and nzn = 0 iff there exists z 6 such that Z = ^ z. 

8 . For a € set (jj{a)b := axb—the cross-product by a on the left. Then a/\b = uj{axb) 
and ALo{a)A~^ = uj{Aa). 

9. If n = 2, then a /\b = b - {Ja)J where J is counterclockwise rotation by 7r/2. 

Proof. All properties are proved by straightforward calculations. □ 

3.2 Left-invariant Riemannian metrics on SE{n) 

Let (•,•) be a left-invariant Riemannian metric on the connected Lie group G with Lie 
algebra g. Let V be the associated Levi-Civita connection. Define R : g x g ^ g by 

(B{u,v),w) = ([u,^],?/). 

If X, Y, Z are left-invariant vector fields on G such that Xg = u,Ye = v, Zg = ru, then from 

(3.1) nVxY,Z) = -{[Y,ZIX) - ([V,Z],y) + {[X,YIZ) 
we obtain 

(3.2) {YxY)e= ]^{[u,v]- B{u,v) - B{v,u)} . 

A left-invariant vector field A is a geodesic vector field if and only if 0 = VjcAT = -B{X, A). 
It is not difficult to show that if the metric is bi-invariant then B{u,u) = 0 for all u € g. 
We adopt the notation: If u e TgG, then g~^v := (dLg-i)^ v 6 TgG = g. 

Proposition 3.2. Let g{f) be any smooth curve in G and X a vector field along g{t), 
not necessarily left-invariant. Define z{t) '■= g{t)~^g{t) and w{f)’■= g{t)~'^Xg(j^). Then 

w + ^ ([^, w] - B{z, w) - B{w, z)). 

\ dt ) g{t) 2 

In particular, g{t) is a geodesic if and only if z = B{z,z). 
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Proof. Let ei,..., e„ be a basis of g and Ei,..., En the respective left-invariant vector 
fields on G. We write X = T.fjEj, g'{t) = Y,hj{t)Ej{g{t)). Then, using Equation 3.2 




Ejfk ^ 2^^ ([L'j) -Lfe] - B{Ej,Ek) - B{Ek,Ej)) 


ait) 


from which we obtain the desired expression after left-multiplying by g{t) □ 

Let the Lie algebra of G be g = s ffi r, where c is an ideal and s is a Lie subalgebra. Let 
(•,•) be a left-invariant Riemannian metric on G and V the corresponding Levi-Civita 
connection. We suppose that s and r are orthogonal subspaces. 

Proposition 3.3. The following properties hold, where we indicate by the same letter 
elements of g and the associated left-invariant vector fields on G. For z,Zj e r, Z,Zj es 

1 . SI z^Zi and B(fZ\,Z< 2 f) Heins 

2 . SJ ziZ 2 e Ifis abelian, V ziZ 2 = 0 and B{zi,Z 2 ) es. 

3. SIzZ = 0 and Szz = [Z,z^. Moreover B{z, Z) s. t and B{Z, z) = 0. 

Proof. All properties follow from the definition of B, Expression |3. 1| for the Levi-Civita 
connection, and the assumption that the subalgebra s and the ideal r are orthogonal. □ 

Let S and R be the subgroups of G having lie algebras s and r, respectively. Then G is 
the semi-direct product G = S x R, where i? is a normal subgroup of G. We now assume 
that i? is a vector subgroup, hence abelian, and that S' is a compact subgroup acting on 
R by linear transformations, S c GL{R), preserving an inner product (•,•) on R. That 
is, S is a subgroup of the orthogonal group 0{R, (•,•)). Elements of G will be denoted 
(A, a) where Ae S and a e R. Indicating the action of S on i? by Aa, the multiplication 
in G takes the form 

(Ai, ai)(A 2 ,02) = (A1A2, Aia 2 + oi). 

Note that Ad(^ o)(0,2:) = (0,Az) and ad (2 Q)(0,2;) = (0, Zz), where {Zz,z) = 0 since A 
acts on R by isometries. 

Proposition 3.4. Under the just stated assumptions B{z,z) = 0 and B(z,Z) = -Zz. If 
g{t) is a geodesic, writing {Z{t), z{t)) = g{t)~^g{t) we have Z = B{Z,Z) and z = -Zz 

Proof. These simple remarks are consequences of the definition of B, the algebraic 
assumptions about the group, and Proposition |3.2| □ 


4 Newtonian mechanics of rigid bodies 

We give here some alternative expressions of Newton’s equation of motion. The approach, 
if not the notations, is essentially that of [1]. Other useful references are [2] and [^. 
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4.1 Momentum of a tangent vector and the momentum map 


If M is a Riemannian manifold with metric (•, •) and v e TgM, we denote 

y{q,v):={v,-)g€T;M 

and call this covector the momentum associated to the (velocity) vector v at (configuration) 
q. The pair {q,v) will be called a state of the system. We often indicate states by 
{q,q), dotting the quantities of which time derivative is taken. The momentum map 


: TM ^ g* was defined in Section 2.3 


If M = G is endowed with a left-invariant Riemannian metric, the momentum map 
for the left-action of G is given by y^{g,g){u) = {g, {dRg)^u)g. Because the Riemannian 
metric is left-invariant, 

3^®(5>5)(m) = (5“^g,Adg-iu)e. 

In this case we also write CP®(u) := (u, •)e e 0 * for u € 0 . Then 

y^{g,g) = Adl^,y^{g-^g) 

where Ad*a = a o Ad^ and Adg is the differential of the map Lg o Rg-i. 

Proposition 4.1. Given a smooth curve g{t) in G and setting z{t) := g{t)~^g{t), then 

jy^{g,g) = Adl,y^{z-B{z,z)). 

In particular, g(t) is a geodesic if and only if momentum y^{g,g) is constant. 

Proof. First note that 

—Adg-i(M) = -[z, Adg-iu]. 

It follows from the definitions that 

= ^(z,Adg-iu) = (i,Adg-izi)- (z, [z,Adg-i]) = {z-B{z,z),Adg-iu). 

The expression on the far right is now lP®(i - B{z, z)) o Adg-i evaluated at u. □ 

When G = SO{n), define on M(n,K) the bilinear form 

(A,y)o:=Tr(Wrt). 


Then (•,-)o is an Adc-invariant non-degenerate positive bilinear form on so(n) and the 
associated left-invariant Riemannian metric on G is bi-invariant. Thus for any left- 
invariant Riemannian metric (•,•) there must exist a linear map £ ; 0 ->• 0 , symmetric 
and positive definite with respect to (•,-)o, such that (^ 1 ,^ 2 ) = |('C(Zi),^ 2 ) 0 - We are 
interested in such L that arises from a symmetric matrix L e M(n,K) according to the 
definition II{Z) := ZL + LZ, in which case 

h:v{L{z,)zl) = yv{z^Lzl). 

If Ml, ..., is a basis of R” of eigenvectors of L, Lui = Aui, then the Ui a Uj comprise a 
basis of so(n) such that L{ui a Uj) = (Ai + Xj)ui A Uj. 
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Proposition 4.2. Given L e M(n, M) define the linear map L : M(n, K) ^ M(n,M) by 
£j{Z) = ZL + LZ. If L is symmetric and non-negative definite of rank at least n - 1, 
then L is an isomorphism and {Zi^Zf) ■= |Tr(£(Zi)Z|) is a left-invariant Riemannian 
metric on SO{n). The tensor B for this metric is 

B{ZuZ2) = [LZi,Z2 ]L-^ 

for all Zi, ^2 6 so(n). 

Proof. Let Ai,..., A; be the distinct eigenvaiues and Vi,... ,Vi the respective eigenspaces 
of L. Let TTj : K" ^ Vj denote the orthogonai projections. Then iTjL = Lttj = XjTTj. It 
suffices to show that L has trivial kernel. Thus suppose Ai(Z) = 0. Then for all i,j, 

0 — TT iTj ^ Z j j — 2 1/ Z TTj 't" 7r 2 Z Ij'IT j — ^ A 2 A J j Hi Zi IT j . 

But Xi + Aj > 0 by the assumptions on L so all blocks HiZiTj are zero, hence Z = 0. The 
expression for B follows from Tr ([^ 2 , = Tr (^[LZi, Z 2 ]L~^Zl). □ 


4.2 Kinetic energy metrics on SE{n) for rigid bodies in 

The left-invariant metrics on G = SE{n) of interest here are derived from mass distri¬ 
butions on the rigid body. Let B c M" denote the body in its reference configuration. 
The position of material point 6 € i? in the configuration g = (^4, a) ^ G = SO{n) x K" 
is $(g,6) := Ab + a. We call ^ : G x B R” the position map and use the alternative 
notations ^{g,b) = g{b) = ^big) as convenience dictates. For now (until we consider 
collisions shortly) B may be any measurable set with a finite (positive) measure p defining 
its mass distribution. Recall from Section 2.3 that m := pi{B) is the mass of the body and 
the first moment of p is 0. When considering the motion of several bodies, we assume 
that the center of mass of each of them in the standard configuration is at 0. 

Elements of g = so(n) x R" will be written in the form ^ = {Z, z). Let Lg and Rg denote 
left and right-multiplication by g. We will very often use the identification G x g = TG 


given by {g,f) 1 -^ {dLg)e£,. Each v € TgG gives rise to the map I 4 • B 


defined by 


Vv{b) = {d^b)qV, which is the velocity of b in state {g,v). The kinetic energy Riemannian 
metric on G is defined so that the inner product of u, u 6 TgG is given by 


{U,v)g= [ Vu{b)-V^{b)dp{b). 
J B 


Proposition 4.3. The Riemannian metric on SE{n) associated to the mass distribution 
pL is invariant under left-translations. 

Proof. To see this, note first that 


K(6) = 

as 


ge"^b= — 
s=o ds 


(Ae^^b + A 

f e*^zdt + a] 

s=0 V 

Jo J 


Here we have used the form of the exponentiation in SE{n) given in Section 3.1 Therefore, 
as A leaves invariant the standard inner product in R", 


Vu{b) ■ K(&) = {Z^b + 0“) • {Z^b + z^) 

and so {{dLg)e^,{dLg)eg)g = {^,f])e- 


□ 
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Recall the inertia matrix L introduced in Section 12.31 
Proposition 4.4. The matrix L associated to mass distribution fj, satisfies: 

1. For arbitrary nx n matrices Zi and Zi, /^{Zib) • {Z 2 b) dfi{b) = Tr . 

2. If A e SO(n), the inertia matrix of the rotated body gB is := ALA'i'. 

3. L = XI if g. is SO{n)-invariant. If g is uniform on a ball of radius R, X= 

Proof. These are obtained by elementary calculations. □ 

Let £j{Z) = LZ + ZL where, from now on, L is an inertia matrix. 

Corollary 4.1. The kinetic energy Riemannian metric can be written in the form 


(4.1) 


{u,v)g = m 


-Tv[lL{Zu)Zl) + Zu-Zy 


where u,v € TgG and their left-translates to g are indicated by and (Zy,Zy). 

Proposition 4.5 (Tensor B for se(n)). Let SE(n) be given the left-invariant Riemannian 
metric associated to the inertia matrix L. Then 

Zi), {Z2 t Z 2 )) = ^^[LZi, Z 2 ] - -Zi A 22^ T ^,-Z2Zi'j . 

Proof. Observe that 


([(Z2,Z2), (^3, ^3)], (^1,^1)) = (([^2, ^3], ^ 2^3 -^3^2), (^1,^1)) 

= m {Tt {[Z 2 , Z 3 ]LZI) + {Z 2 Z 3 - Z 3 Z 2 ) ■ zi} 

= mTr Z2] - -zi a z 2 j j - m(Z2Zi) • Z3 

= -^ 2 ] - - 2:1 A Z2 j T -Z 2 Z 1 j , {Z 3 , Z3)| . 

The claimed identity now follows from the definition on R. □ 

We note that if L = A/, then B{{Z,z),{Z,z)) = {0,-Zz). 

Proposition 4.6. Give G = SEigT) the left-invariant Riemannian metric defined by a 
mass distribution on the rigid body B with inertia matrix L and mass m. Let f=(lT,'u;)€g 
and (g, v) € TG where g = (A, a) and v = {dLg)e{Z, z). Then 

5’®(ff,w)(C) = [{MaE{Z) + Xc/\Vc)W^] + mvc-w. 

Here Xc = a is the position of the center of mass of the body in configuration g and Vc ■= Az 
is the velocity of the center of mass for the given state {g,v). 

Proof. This is a straightforward computation based on the definition of CP® and the 
expression of the Riemannian metric given in Corollary | 4 . 1 [ □ 
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Let (•,-)g be the left-invariant inner product on g given by 


(4.2) ((Z, w) )g := Tt{ZW^) +z-w. 

Then, with the notation of Proposition |4.6[ 

(4.3) (AdA£(^) +a:c AUc) ,Ucj , {W,w)^ 

4.3 Singular force fields and impulses 

Let M be the configuration manifold of a mechanical system with the kinetic energy 
Riemannian metric and material body B with mass distribution measure /i. 

A force field on M is a bundle map F -.TM ^ T* M possibly depending on time, although 
we omit explicit reference to the time variable. So if e T^M, then F{q,v) e T*M. Forces 
acting on B typically arise from a K^-valued (possibly time dependent) measure ip on B, 
parametrized by TM , called the force distribution. From such a measure we define the 
force field F{q,v) eT*M such that for each ueTqM, 

F{q,v)(u)= f Vu{b)-dpq^y{b). 

•J B 

We are interested in cases where is singular, supported on a single point of B. 

Definition 4.1 (Newton’s equation). Newton’s equation of motion of the (unconstrained) 
mechanical system with configuration manifold (M, (•,•)) and force field F is 

^^(9,9) = F{t,q,q). 

Proposition 4.7. Give M = SE{n) a left-invariant Riemannian metric and let F be a 
force field on M. Let F"^ be the dual field, so F{t,q,v){u) = {F'^{t,q,v) ,u)q. Then 

^w-B{w,w) = {dLg-i)gF* ^ jy^{g,g) = R*gF. 

Proof. These are consequences of Propositions |3.2| and |4.1[ □ 

When M = SE{n), it is useful to regard the force field as a Lie algebra-valued by 
left-translating each force vector to TeG. We define 

‘^{t,g,g) := {Y{t, g, g),y{t, g, g)) := {dLg)l^F*{t,g,g) eg. 

Then, using the notation of |4.2| 

(T,Adg-i^)g = ]^’Tr[{NdAL{Y) + Xc A Ay)W'^] + {Ay)-w 
( 4 - 4 ) ^ , 

= (AdA£(y) +XcAAy),Ayj,{W,w)j 
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Proposition 2.2 follows from these remarks, keeping in mind that Vc = Az and Xc = a. 

Let the force be applied on a single point Q = Q{t,g,g) € M", so that (p is supported on 
Q. Let u = (dRg)^ 6 TqM, ^ = {W,w) 6 g. Note that we are using right-translation here 
since we wish to evaluate the last of the equivalent equations of Proposition |4.7| on the 
Lie algebra element Then there exists 3 = 3(t,g,g) depending only of F such that 

F{t,g,g){u) = f Vu{b) ■ dpg^g{b) = Vu{g^^Q)-3= {WQ + w)-3= /\Q)W‘^) + 3-w. 

J B A 

This gives 

It follows from the expression |4.4| of F that 


Ady 4 £(y) + Xc A Ay = Q A J 
Ay = 3. 


Writing /c := Ay, this is equivalent to fc = 3 and Adyi£(y) = {Q- Xc) a /c. Therefore, the 
equation of motion becomes 


mvc = fc 

mAdA {d3Z - [LZ, Z]) = {Q - Xc) A fc 


proving Proposition 2.3 


Our informal discussion of the idea of impulse from earlier in the paper and the above 
remarks now give the expression ((Q - Xc) A 3c, 3c) 6 g for the change in momentum due 
to singular forces applied to Q. This gives the following. 


Proposition 4.8 (Change in momentum due to impulsive forces). If the rigid body with 
mass m, inertia matrix L and associated Lie algebra map L, is subject to an impulsive force 
concentrated at point Q € M" at a given time, then momentum changes discontinuously 
according to 


(4.5) 


mv)i - mv^ = 3c 

m Ady 4 C(Z'^ - Z~) = {Q- Xc) A 3c 


for some vector 3c 6 R" depending on the state of the body. As before, Xc = a indicates 
the center of mass of the body in configuration g = {A, a), = Az"*" are the velocities of 

the center of mass immediately prior to and after the application of the impulse, and 
{g, (Z"^, z"^)) are the pre- and post-impulse states of the body. 


4.4 Several bodies and momentum conservation 

If the mechanical system consists of several unconstrained rigid bodies, Bi,... ,Bk (in 
reference configuration) subject to forces Fj{q,v), j = I, ... k, the configuration manifold 
M is a subset of the product G x ••• x G, with one copy of G = SE{n) for each body. We 
consider for now only motion in the interior of M. 
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We say that forces are internal to the system if they are somehow due to the influence 
of the bodies on each other. More specifically, we use term ‘internal’ when Fj = Fij 
and the F^j = Fij{q,v) —the force body i exerts on body j in state {q,v) —satisfies the 
property of action-reaction: Fij = -Fji. If the forces are derived from, possibly singular, 
measures ipq^v{j,b\i,b') on Bi x Bj so that 

Fij(q,v){u)= f f Vu(b')-d:pq^y{j,b'\i,b), 

then the action-reaction property, expressed in terms of </?, means that 

d‘i:>{i,b\j,b') = -dip{j,b'\i,b) 

for almost every b,b' (with respect to (p). Newton’s equation applied to body j is then 

itj 


and the total momentum CP(g, q) = is conserved: 

9 ) = E E d, q) = 0 . 

at j ^^j 

Another way to interpret the notion of forces internal to the system is to assume that 
the total work the Fi do along a rigid motion of the entire system, that is, the work along 
a path in M of the form 7 (t) := = (e^^qi,..., e^^gu) is zero. The total work is then 

0= r Y.^3ili^)n'{t)){ij{t))dt 

3 

= E f\R;F,i^{t),l'{tmOdt 

3 

= E/ J^[‘^jil3(^)^^j(t))(0]dt 

= E5’®(^W’V(fe))(0-E5’”(7(a),7'(a))(0 

3 3 


and, again, the total momentum (now in the sense of the momentum map on g) is constant. 
In this sense, conservation of momentum follows, as expected, from a symmetry property. 

Of particular interest here are two bodies that interact through impulses applied to a 
common point of collision Q. Then for each body, indicated by the index z = 1,2, 


(4.6) 


ruivF - niiV^ i = 3c,i 
TTli A.(ij^^dLi(^Zi — Zi ) — (^5 ~ 3Cc,i) 3\ 3c,j 


where the impulse vectors satisfy 3c ,1 + Uc ,2 = 0 by conservation of momentum. Proposition 
|2.4| is now a consequence of this observation and of Proposition |4.8[ 
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5 Kinematics of two rigid bodies 


The configuration manifold of a system of several (unconstrained) rigid bodies in M" is 
a submanifold with boundary of the product of copies of the Euclidean group SE{n), 
one copy for each body. The Riemannian metric is then the product of the Riemannian 
metrics for each single body. Here we focus on the boundary of the configuration manifold 
of two bodies and certain structures therein. 

Let Bi and B 2 be submanifolds of M" of dimension n having smooth boundary and 
equipped with mass distribution measures /Zi and ^2 with masses rrij := fj,j{Bj) < 00 and 
zero hrst moment. The bodies need not be bounded. The configuration manifold M is by 
definition the closure of 

(5.1) Mg := {<7 = (gi, ( 72 ) € G X G '■ gi{Bi) n g 2 {B 2 ) = 0 } 

where G = SE{n). We further assume that each collision configuration q = ( 771 , 772 ) e dM 
is such that 771 (i?i) n 772 (^ 2 ) consists of a single point. 

The definition of M as the closure of Mg is not a very useful description of M near its 
boundary. In particular, it is not so clear how to translate geometric information about 
the boundaries of the Bj into information about the boundary of M. For this purpose we 
introduce the extended configuration manifold Mg defined below. 

Let Nj be the boundary of Bj and let Vj be the outward-pointing unit normal vector 
field on Nj. By a (positive) adapted orthonormal frame at 5 6 Nj of sign e 6 {+,-} we 
mean a positive orthogonal map a ; M" ^ TfelR" = K" such that (Te„ = evj{b). Here e„ is 
the last vector of the standard basis (ci,..., e„) of M”. Hence a is an element of SO{n) 
mapping M" isometrically to T^Nj. If a is an adapted frame and h e E[ := SO{n- 1), then 
ah is also an adapted frame with the same base point as a. In this way, E[ acts freely and 
transitively by right multiplication on the set of adapted frames at any given point of Nj. 



Figure 13: Interpretation of the map 'P. The transformation sends body Bj from its standard 
configuration to the configuration that takes the adapted frame Uj to the standard 
frame in R", and the point bj into the line through the origin along e„ a distance 
s/2 from the origin. 


We denote by Tj the principal iL-bundle of adapted (positive) orthonormal frames over 
Nj of sign e. Elements of T) will be written {b,a), where cr is in the fiber 3^j{b). The 
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extended eonfiguration manifold is the product Mg •= x 5’2 x G x [ 0 , oo). We can now 
define the map : Mg -> G x G by 'i>{bi,< 7 i,b 2 ,< 72 ,g,s) = (31,32) where 

3 1 = 33 i = - s/ 2 ) = |Ae„ - 

32 = 332 = g { f ^ 2 ^,-(^ 2 ^ b 2 + sl 2 ) = ^ Aa 2 ^, a + 

The geometric interpretation of T is shown in Figurej^ Note that points on the boundary 
of M correspond under T to points in Mg with coordinate s = 0 . The groups G and H 
naturally act on Mg on left and right, respectively: 

3(61,0-1,62,0-2,3': •= {bi,crih,b2,(T2h,ggh,s). 

The quotient Mg/M is easily seen to be a smooth manifold and the projection Mg ^ Mg/M 
is a principal M-bundle. It is also immediate from the definitions that 

'i’igqh) = 3^(3) 

for all f 6 Mg, where the action of G on GxG is defined by 3(31,32) = (331,332)- Therefore, 
induces a G-equivariant map 

T : Mg/M ^ G X G. 

Equivariance means 'I'(33) = 3'!'(3). The G-action on Mg admits a smooth cross-section: 

S ■■= 3^1 X J2 X {e} X [0, 00). 

The G-action on Mg and on Mg/M leaves invariant the coordinate s; in particular, it 
leaves the boundary of these two manifolds invariant. 



Figure 14: For M^jH to be a good parametrization of M near the boundary some pathologies 
must be avoided. Far left: a boundary configuration in Me/M that is not in SM; 
middle pair: two distinct elements of Me/M corresponding to the same element in 
SM; far right: a curve in Me that is mapped under to a single point in M. 


It is natural to expect that under reasonable assumptions the restriction of di to a 
neighborhood of the boundary of Mg/M will be a diffeomorphism onto a neighborhood 
of the boundary of M, thus providing a useful parametrization for the purpose of 
understanding collisions. Figure [T 4 | shows some of the situations that must be avoided. 
We give shortly a few sets of sufficient conditions for T to be a local diffeomorphism, but 
our immediate goal is to explore Mg, Mg/M, and their boundaries a little further. 
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5.1 The tangent bundle of dMe and d { MelH ) 

Recall that the shape operator of a hypersurface N c K" with unit normal vector field i/ 
at a point b e N is the linear map Sb ■ TbN TbN defined by u i-*- -D^v where Dy is the 
Levi-Civita connection for the standard Euclidean metric in K". We write VyX := libDyX 
for a tangent vector field X on N ^ where Ilf, is the orthogonal projection from K" to TbN . 
This is the Levi-Civita connection on N for the induced metric. Let (6,(7) be a point in 
the adapted frame bundle T over N and {v, Q a tangent vector to at {b, a). Then, by 
differentiating = -(-l)'^cr(t)e„, where is a smooth curve representing 

(u,C) at {b,a) = (7(0),cr(0)), we obtain 

Sb{v) = (-l)Ten = -Co-~T(5) = -aVa~^v{b) 

where V := cr~^C can be regarded as an element of so(n) just as a is viewed as an element 
of SO{n). The tangent bundle of T for any smooth hypersurface N has now the following 
description. Let {b,a) e Then 

T(f, = {(u, V) € TbN X so{n) : Sb{v) = -aVa~^i/{b)} . 

As before, we use the canonical identification TG = G x g and write elements of g in the 
form {Z,z) 6 so{n) x M". The shape operator of Nj will be written We omit the 

superscript when it is clear from the context to which body the operator is associated. 
Then the tangent space of Me at a point q = {bi,ai,b2,cr2,g,s) is given by 

TqMe = {{vi,Vi,V 2 ,V 2 ,Z,z,g) : Vj € Tb.Nj,Vj €S 0 (n), (Z,z) € g, M, 

Sb^Vj = -crj V,cr"Tj(6j),j = 1,2} . 

Tangent spaces to dMe consist of those vectors for which g = 0 . 

Let Gq and Hq represent the orbits through q 6 dMe of the (right and left, respectively) 
actions of G and H on Me- The tangent spaces at q of the respective orbits will be 
written g^ and f),j. Then 

( 5 . 2 ) g, = {( 0 , 0 , 0 , 0 , Z, z, 0 ) : {Z, z) 6 g} and f), = {( 0 , Y,0,V,V,0,0) : Y e t)} . 

At any q e dMe the differential of 'h is 

d^q{vi,Vi,V2,V2,Z,Z,Q) = ( , Zi, ^2, Z2 ) , 

where, denoting Adcr(IT) := crITcr“^, 

Z. = kd^ AZ -VX 

( 5 . 3 ) ^ 

Zj = ajZ - Ado-^ {Z - Vj)bj - vj - gvj{bj). 

The next proposition contains as a special case Proposition | 2 . 1 [ It uses the notation 

Sj := ^ 

for any q = {bi,ai,b2,cr2,g, s). We allow s to be non-zero, in which case is the shape 
operator of the level hypersurface of MejH corresponding to value s. 
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Figure 15: Situations for which Proposition 5.1 applies. 


Proposition 5.1. The map 'E': ->■ G x G is a submersion from a neighborhood IL of 

the boundary of Mg onto a neighborhood of the boundary of M if any of the following 
conditions involving U and the shape operators holds. 

1. Si + S 2 is non-singular at all points of dM^,. In this case, U is a neighborhood of 
dMe where this non-singular condition holds. 

2. U is a neighborhood of dMg where one of the Sj is non-singular and Si+ S 2 - SS 1 S 2 
is also non-singular. 

3. If the two bodies are convex and the boundary of one of them has non-vanishing 
Gauss-Kronecker curvature so that Sj is everywhere non-singular on Nj for some j, 
then li = Me. 

In each case U is G-invariant, the kernel of d'i’q is t)q at each a ill, and 'E'ju • U lE'(ll) 
is a principal H-bundle. In addition, the boundary of M is a smooth submanifold and 
there are smooth functions bj ■ DM -»• Nj, j = 1,2, such that bi{q),b 2 {q) are the unique 
points on the respective bodies that are brought into contact in collision configuration q. 


Proof. By counting dimensions we see that 11 should be a neighborhood of the boundary 
of Me where the kernel of d'ltq is l)q. It follows from equations 5.2 and 5.3 above that 


this kernel contains 1}^. We show equality under the conditions of item (2), the other 
cases being similar. Say that S 2 is non-singular. From the explicit form of d'l'q given 
in|5.3 we see that f = {vi,Vi,V 2 , V 2 , Z, z, g) lies in that kernel if and only if Z = Vi = V 2 
- (-l)^ |Ze„ for j = 1,2. Observe that 


and z = a, 


"Vj and Zcr, 


lie in 


which is 


orthogonal to e^. Hence g = 0. Keeping in mind UjVjCn = {-lySb Vj, we obtain 


-SiZ — ZCji ^SiZCn 

Q 

S2Z = ZCn — —S2Zen. 


From this we conclude that [^i + S 2 - sSiS 2 ]Sf^Zcn = 0 which, under the conditions of 
(2) implies that Ze„ = 0. Since S 2 is non-singular, this also implies that z = 0 and Vj = 0. 
Therefore, f = (0, Z, 0, Z, Z, 0,0), where Z 6 () since Ze„ = 0. That iklu is a principal 
H-bundle is now easy. G-equivariance of ik implies that U is G-invariant. □ 
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Assuming for simplicity that U of Proposition 5.1 is all of Mg (we only need what 
follows on some neighborhood of the boundary of Mg), it is useful to know whether the 
principal bundle Mg ->• M admits a G-invariant connection since the associated horizontal 
subspace ‘Kq can then serve as a proxy for the tangent space of TgM, without having 
to go to the quotient. A principal i/-connection on Mg is given by a one-form w taking 
values in f) and satisfying the properties: 


1 . ujq{Yq) = Y e i), where Yg is the vector induced by the infinitesimal action of (); 

2. h*uj = Adfi-i o uj. 

Let n be the orthogonal projection from R" to = e^. 


Proposition 5 . 2 . For any real constants ci, C2, C3 such that ci + C2 + C3 = 1 the Ij-valued 
one-form uj on Mg given by 


ujq{vi,Vi,v2,V2, z, z, g) = ciiii/iii + C2nv2n + C3nzn 
is a G-invariant H-connection on Mg. 

Proof. This is a simple check. □ 

Let us choose ujq{f,) := nViII and denote by iKg the horizontal subspace defined by this 
choice of connection form. Recall the maps Sj on for q = (61, cti, 62, < 72 , g, 0 ) € dMg 
given by Sj = af^Sb^crj- The hj are determined uniquely from 4 '((?) and the aj are 
determined uniquely up to an overall common element of H acting on the right. 

Proposition 5 . 3 . Let q = (bi,ai,b2,< 72 , 9 ,d) ^ dMg and suppose that Si + S2 is invertible. 
Then d'i’q maps TCg isomorphically onto 0 x g and 


d'dtq (TCq n Tq(dMg)) - {(Zi, zi, Z 2 , 22) 6 0 X 0 • • (^1^1 + zi) + 1/2 ■ (^2^2 + Z2) - 0 } 


where vi = I'libi) and V2 = V2{b2) are the unit normal vectors. 

Proof. The proof is elementary, but we show the main point. Let f = {Zi, zi, Z2, Z2) e g xg. 
We wish to show the existence of a unique f = (vi, Vi, W2) 1^2) Z, z, g) € TCq that is sent to f 
under dg'dt. The components of ^ satisfy: 

afvj 6 ny,n = O, {Z,z) eg, p e R, S^a^Vj = {-lyV^en- 

and Zj , Zj are related to theses quantities by 

Zj = Ado-j {Z — Vj), Zj = ajZ — Ado-^ {Z — Vj)bj - Vj — gvj{hj). 

Writing g = {A,a) and ^'(g) = (51,52), we have gj = {Aa~^,a- Aaf^bj). Note that 
(5.4) crl^ui - cr2^U2 + 2gen = cr2H^2&2 + Z 2 ) - ay{Zibi + Zi) 
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from which we obtain g and ai^vi - a2^V2 ^ R" ^ in terms of the Zj and Zj. Also 

- Ad^-iZi j e„ = {V2 - V2)en 

= SiCJ^^Vi + S 2 (J 2 ^V 2 

= {Si + S2)<Ji^Vi + S2{(J2^V2 - CTl^Vi) 

= (^1 + S2)(T-i\l + S2li-{(T2\Z2b2 + Z2) - (J-^^Zibi + Zi)] 


from which we obtain CTi^ni in terms of the Zj and Zj under the assumption that + S2 
is invertible. From Proposition 3.1 item (8), we deduce 

Vi = nViII + e„ A Vie„ = e„ A ViCn = -e„ a Sia^^vi 


so that Vi is also uniquely determined by the Zj and Zj. From 

P2-lA = Ad^-iZ2-Ad,-iZi 

we obtain V2 uniquely and from the above | 5 . 4 | we obtain V2 uniquely. From these we 
easily obtain Z and z as well, proving the first part of the proposition. The second part 
follows from the observation that a vector is tangent to dM if and only if = 0 . □ 


5.2 The non-slipping, rolling, and diagonal subbundles 

Let 7(t) be a smooth curve in dMg such that q = 7(0) = (&i,(71,62,0-2,g, 0 ). We omit the 
variable s, which is set to 0 for a boundary point. Let ( 7 i(t), 72 (i)) be the image of 7 
under 'll and write 7 ^( 0 ) = gj, q'■= (31,(72), where gj = (Aj,aj) and g = (A,a). Denote 
the components of the infinitesimal motion in Mg by 


^= 7 '( 0 ) = (?;i,Ui,(;2 ,T2 ,^,z), 


omitting 3 = 0. The two bodies in configuration q are in contact at 31(61) = 32(62). 

Definition 5.1 (Non-slipping and non-twisting conditions). The infinitesimal motion 
^ € TqMf, is said to satisfy the non-slipping condition if the velocities of the material points 
bj at the contact configuration are equal. It is said to satisfy the non-twisting condition if 
the tangent planes to Nj at bj do not rotate relative to each other under T 

We now derive an explicit expression for these conditions. The infinitesimal motion of 
Bj is given by £^j := {Zj,Zj) € g, which is obtained from d'^tq^. We know that 


Z,=Ad,^{Z-Vj) 

Zj = ajz- (Z-Vj)bj - Vj. 


Due to Proposition 4.3 V^. {bj) = Aj{Zjbj + Zj) = A(z- cr ■ ^Vj). The non-slipping condition. 


V{i(6i) = ^{2(^2), then reduces to 


( 5 . 5 ) 


(TjTi = a 2 ^V 2 . 
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Turning now to the non-twisting condition, let Uj be a tangent vector to Nj at bj such 
that, in the contact configuration given by q, is sent to a common vector, for j = 1,2, in 
the plane of contact. Thus AiUi = A2U2- The infinitesimal rotation of AjUj at the point 
of contact is 




■ Vj)aj Uj. 


The orthogonal projection to the plane of contact is AIIA~^, recalling that II is the 
orthogonal projection to (It may be helpful to keep in mind Figure 13 ) Because 

A = Ajaj, the non-twisting condition takes the form nViII = nV2n and since nkill = 0 
holds for horizontal vectors, nVjTI = 0 for j = 1,2. 

Now let 'i’iq) = (51,52), 5j = (Aj,aj) and ^ = {Zi, zi, Z2, Z2) = d'Sq^. The non-slipping 
condition expressed in terms of ^ becomes 


Ai[Zibi + 01 ] - A2[Z2b2 + Z2^ 
and the non-twisting condition becomes 


AdAjZj = W + i^j{bj) A Wj 


for a kF 6 so{n) independent of h and Wj € T^-Nj. 


Definition 5.2 (Non-slipping, rolling, and diagonal subbundle). The non-slipping sub¬ 
bundle ofT{dM) consists of all tangent vectors satisfying the non-slipping condition. The 
rolling subbundle ofT{dM) consists of all tangent vectors satisfying both the non-slipping 
and non-twisting conditions. The diagonal subbundle ofT{dM) is the tangent bundle to 
the orbits of the action of G on DM defined by 5(51,52) = (551,552)- We denote these 
three subbundles, respectively, ©, 91 ,S). We refer to these collectively as the kinematic 
subbundles of T{dM). Notice that = Qq, using previous notation. 


Starting from this definition rather than Definition 2.1 the content of the latter becomes 
a statement, which is proved by the above remarks. 


6 Collision maps 


Let now M c G x G be the configuration manifold of two rigid bodies in K", where 


G = SE{n). By condition 2 of Proposition 5.1 M has smooth boundary and boundary 
points represent configurations in which the bodies are in contact at a single point. Let 
the state of the bodies before and after collision be given by the element of TqM, q € dM, 
represented by 

(Zf , zt , Z^ , Z2) egxg. 


Here the sign ‘+’ indicates post-collision velocities and pre-collision velocities. We 
obtained in | 4 . 6 | a condition on the pre- and post-collision velocities due to impulsive forces 
that act at a single point of the body. We restate it here. Let the common point of 
contact he Q = Ajbj + Uj, where bj is the material point in standard body configuration 
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which corresponds to Q in configuration gj = where aj = Xcj is the center of 

mass of body j in the given configuration. Then we obtain from expression |4.6[ 


( 6 . 1 ) 


4 

Zj = Zj + ^ {bj A Uj ) 


where Uj = Aj^dcjlmj. We should add to these equations Jci + Jc 2 = 0 for conservation of 
(linear) momentum. 


Proof of Theorem 2.1 A simple dimension count gives dim (Lq = n and dim &q = 2 dimg-n 


so that the sum of the two dimensions equals dim T^M. Therefore, it suffices to show 
that these subspaces are orthogonal. The Riemannian metric on M is the restriction of 
the product metric on G x G (each factor having a possibly different metric as the bodies 
may have different mass distributions.) Explicitly, let u,v € TqM and write 


= {{Yi,yi),{Y 2 ,y 2 )), w = ((Zi,zi),(Z 2 ,za)). 


Then 

( 6 . 2 ) 


^)9 = E’ 

3 




Now consider the vectors 

V = A Ui),Mi), {Ti2^{b2 /^U2),U2)) € 

w = {{A\^ZiAi,Al^z* - Al^ZiAibi),{A-2^Z2A2,A2^z* - A2^Z2A2b2)) e 6 , 
where Zj = Z - Vj and z* = z - z'. Observe that 

Tr ((&j A Uj)(Aj ZjAj)"^)^ = 2uj • {Aj ZjAjbj). 

Then 


{v,w)q = -Tr((&j auj)(A/ZjAj)1')) + {A^^z* - A^^ZjAjb^) ■ uj 

3 ^ 

^ X! ['^3 ’ (^j ^3 ^3^3 ) ■*■ ^ ^3 ^3^3 ) ’ '^j] 

3 

= J^mj(Afz'^)-U3 


3 

= Z 


Y^mjAjUj. 

j 


But miAiUi + 7712^2142 = 0 by the definition of <Lq so the two vectors are orthogonal. □ 


Proof of Corollary 2.3. Let G be a linear involution in O(n-l). Then G is diagonalizable 
over K with eigenspace decomposition = (G + ffi (G - and eigenvalues 

1,-1 having multiplicities n - fc - 1 and k, respectively, where k € {0,1,..., n - 1}. Thus 
for each such G there is k and A e GL{n - 1,IR) such that G = A~^JkA where Jk is the 
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diagonal matrix diag(/„_fc_i,-/fc) and Ii indicating the I x I identity matrix. We can 
take A to be orthogonal. In fact , let A = SU be the polar decomposition of A into a 
positive symmetric part S = \/A^ A and orthogonal part U. The condition C^C = I implies 
that and commute, from which it follows that 5'^, hence S, is also a block matrix 
with 0 on the off-diagonal blocks of size fc x (n - fc - 1) and {n - k - 1) x k. Therefore, 
S commutes with Jk whence the claim. Thus the set of all orthogonal involutions in 
dimension n - 1 is the disjoint union of the sets 2k = {U^JkU : U € 0{n - 1)}. It is clear 
from this description that 2k is the homogeneous space 0(n-l)/L, where L is the isotropy 
group of Jk- Equivalently, L is the subgroup of all U that commute with J^, which is 
easily seen to be the product 0{n - fc - 1) x 0{k). □ 

The following proposition gives a concrete expression for the unit normal vector field. 

Proposition 6.1. Let k'j{bj) denote the unit outward pointing normal vector to body Bj 
at the boundary point bj. Then the unit normal vector to dAI at q is given by 

n, = (Ci(£j;^(6i A Vi{bi)),Vi{bi)),C2{L2^{b2 a l^2(&2)), i^2(&2))) 


where Ci,C 2 are defined up to a common sign by the equations miCi = m 2 C 2 and 


Y,c]mj 

3 


1 + ^ A Vj{bj))y) 


= 1 . 


Proof. The unit normal vector n^, being an element of can be written as 

AMi),Ui), {L2^{b2 A U2),U2)) 

for some uj € K". Recall that a vector v = ((Zi, zi), (Z 2 , Z 2 )) tangent to DM has the form 

z,=Ad^yz-v,) 

Zj =(7jZ- (Z - Vj)bj - Vj 

where Vj and Vj are related through the shape operators as discussed earlier and Vj is 
tangent to the boundary of body Bj at bj. Let as before Vjfibj) denote the unit normal 
vector to body Bj at bj. Using the explicit form of the Riemannian metric we obtain 
after straightforward computation that 

d={Ta.q,v)q = -Y,mjVj-Uj. 

3 

This being true for all Vj implies that Uj = CjVj{bj). But + m 20 'yu 2 = 0 by the 

definition of and cryvjfibj) = -{-Ifien. Thus the first equation. The second equation 
corresponds to the condition ||ngp = 1. □ 


7 Proof of Theorem 2.2 


Define the one-form 9 on TM from the kinetic energy Riemannian metric on M so that 
= {vjdiTy^) for each ^ e Tq^y(TM), where dir^ is the map induced on the tangent 
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space at {q,v) of the base-point projection map tt : TM ^ M. We briefly recall the 
definition of the vertical and horizontal subbundles E'' and of T{TM). For simplicity 
of notation we denote points in TM by v rather than {q,v). Then the fiber E^ above 
V is the kernel of diTy and E^ is the kernel of the connection map Ky : Ty(TM) ->• TqM, 
defined as follows: if ^ = w'(0) where w{t) is a curve through v representing then 
~ ^ ^ vector fields on TM, then 

( 7 . 1 ) dey{X,Y) = {KyX,dTTyY) - {KyY,dTTyX). 

See 0 for more details. 

Now let S denote the boundary of M, N the pull-back to S of the tangent bundle TM 
under the inclusion map and for each value £ > 0 define 

:={(g,^;)6fV:i||u||2 = £}. 

So is a level set of the kinetic energy function. It is shown in |S], that the pull-back 
of dO to under the inclusion map is non-degenerate on \TS, and so it defines 
there a symplectic form. If the ambient space of the system is K” then has dimension 
2n- 2. The canonical billiard measure is now the measure associated to the (2n- 2)-form 
n = pulled-back to \ TS. 

The smooth field q Qq oi collision maps defines a smooth map (away from singularities) 
C : . The pull-back of 9 under this map is easily shown to be 

{e* 9 )y {0 = {v,eqd 7 ry 0 q. 

Note that diTy^ € TqS whenever ^ € TyN. Define the projections 11^ from TqS to the 
eigenspaces of Cq associated to eigenvalues ±1. The assumption that 6 is parallel is 
equivalent to one of these projections (equivalently, both) being parallel. Now define 
9y{^) = {v, IlgdiTy^), so that 0 = + 9~ and Q*0 = 9'^ - 9~. Consequently, 

e*d9 = d9^ -d9-. 

The projections 11'*' can also be defined on TN^ by requiring 

n*d7r„ = d7r„n^, KyIly=IlqKy. 

Using these maps we define 2-forms w'*' by ujy{^,r]) := d9y{]Y^^,Yl^r]). We now wish to 
relate w'*' and d 9 ^ . 

First define a tensor field on S such that for u,v e TqS and any vector fields X, Y 
on S such that Xq = u and Yq = v, we have 

^?*(u,^;) :=(V„n*)U-(V„n^)X. 

It is not difficult to verify that this is indeed a tensor field and the definition does not 
depend on the extensions X,Y of u,v. Furthermore, 'i?'*' vanishes under the conditions of 
Theorem |2.2[ A straightforward calculation now shows that 

d9ti^,v) = w*(C,77) + {v,d^{dTry^,dTryr])). 
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Therefore, dO^ = a;'*' when the field of collision maps is parallel. Moreover, d6 = oj^ + oj 
and 6 * 0 ;''' = ± 0 ;'*'. It is now easy to check that 

= ±(0;*)"" A (0;-)"- 

where n± are the dimensions of the eigenspaces of Qq associated to eigenvalue ± 1 , and 
we finally obtain 6 *(d 6 *)”“^ = ±{d9)'^~^. Therefore, the measure induced by is 

invariant under 6 . 
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